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Abstract 


across the entire width of the flow, is 


A method is presented for estimating the dis- 
persion and dilution of jet aircraft exhaust from 
aircraft passage through times on the order of 
weeks thereafter. In the near wake of the aircraft, 
the solution is that for round turbulent jets in a 
parallel flow. More rapid dispersion due to atmos- 
pheric effects begins when the scale-dependent eddy 
viscosity becomes larger than the turbulent jet 
eddy viscosity. In the far wake region, the solu- 
tion approaches that for scale-dependent dispersion 
from a point source moving with the aircraft. Cal- 
culations are presented for supersonic aircraft at 
high-altitude flight conditions. 

Introduction 
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where Ux is the local streamwise velocity, u„ is 
the flight speed, Uj is the jet exhaust velocity, 
rj is the radius of the jet at the engine exhaust 
plane, and p is the gas density. All symbols are 
defined in the appendix. 

A second relation is obtained by evaluating 
the integral momentum relation at the half -value 
radius, r 5 . 


Concern about the possible adverse Impact of a 
fleet of high altitude aircraft on the environment 
has motivated Interest in aircraft wakes and 
stratospheric chemical kinetics. The significance 
of chemical reactions in aircraft wakes is expected 
to be dependent on atmospheric conditions and on 
the concentrations of possible reactants. The 
latter, of course, is determined by the engine 
exhaust emissions and their dispersion into the 
atmosphere. 

Aspects of the dispersion in regions dominated 
by the aircraft initial conditions have been con- 
sidered in References 1 and 2. The dispersion in 
regions dominated by ambient atmospheric conditions 
has been discussed in Reference 3. The method pre- 
sented in this report extends the analysis of Ref- 
erence 2 to Include atmospheric dispersion effects 
and multiple engine aircraft configurations. At 
long times after aircraft passage, the solution 
approaches the result for scale-dependent disper- 
sion from a point source moving with the aircraft. 

Although the solution approached in the far 
wake region using the present method is similar to 
that given in Reference 3, the present method does 
not require knowing the time after aircraft passage 
at which atmospheric effects become Important or 
the plume size and concentrations at that time. 
Since calculations are begun at the engine exhaust 
plane, the initial conditions required for the 
method of Reference 3 are included in the results 
obtained using the present method. 

Analysis 
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where rj is the radius at which the local ve- 
locity is one-half the sum of the centerline ve- 
locity, Ug, and the flight speed, u®. A gradient 
diffusion model for the shear stress is assumed. 

The eddy viscosity, e, in the near wake region is 
assumed to be directly proportional to the jet mass 
flow plus the absolute value of the entrainment 
mass flow and inversely proportional to the jet 
width. The resultant eddy viscosity relation, from 
Reference 2, is 



The equations are closed by assuming that the 
velocity difference and concentration difference 
ratios are Gaussian, that is. 



The dispersion of the exhaust of a turbojet 
engine at flight conditions is modeled as a round 
jet mixing with a parallel flow. The Prandtl and 
Schmidt numbers throughout the flow are assumed to 
be equal to 1, thus the conservation equations for 
momentum, total enthalpy, and species concentration 
are identical. This permits the problem to be 
solved by considering only the continuity and mo- 
mentum equations. 

The Integral momentum equation, evaluated 


where 1 = In 2. Equations (1) to (4) plus the 
perfect gas law provide a complete set of equations 
for the solution of the jet spreading and dilution. 
The integrals in Equations (1) to (3) may be evalu- 
ated in closed form at any streamwise location. 

The jet width is obtained in terms of the center- 
line velocity with Equation (1) . The incremental 
changes in streamwise distance corresponding to 
selected Incremental changes in the centerline ve- 
locity are obtained with Equation (2) . The stream- 
wise distance, x, corresponding to any desired u^ 



must then be obtained by numerical Integration. If 
desired , the spreading and dilution results may be 
expressed In terms of the time after aircraft pas- 
sage rather than the downstream distance. In this 
case, t = x/u„,, where u,, Is the flight speed. 

At downstream distances where p = p„ and 
(u^ - Uoo) « u„. Equations (1) and (4) provide the 
following relation between the half-value width and 
the centerline concentration difference ratio. 



where c^ is the concentration on the centerline 
at any downstream location. In this region the 
asymptotic result for a small excess jet Is ob- 
tained , namely 
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and the eddy viscosity Is decreasing with down- 
stream distance, l.e.. 


some distance Its value becomes less than the eddy 
viscosity for scale-dependent diffusion in the at- 
mosphere. When this occurs, the atmospheric ef- 
fects are assumed to be dominant and the scale- 
dependent eddy viscosity Is used in all subsequent 
calculations. Up to this point, Che exhaust of 
each engine Is assumed to decay axlsymmetrlcally. 
However, when atmospheric effects are dominant, the 
plume dispersion in the horizontal plane Is known 
to be greater than the vertical dispersion. (3 ,4) 
Accordingly, a blvariant Gaussian distribution is 
assumed for the dispersion of the exhaust from each 
engine. That Is, 



where y^ and Z 3 are the horizontal and vertical 
half -value widths, respectively. For a twin-jet 
aircraft, the concentration distribution In the 
plume is given by 
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Results for the dispersion of both large and small 
turbojet engines from the engine tailpipe through 
the above region are given in Reference 2. 

Multiple engine aircraft are treated by super- 
imposing the concentration distributions from the 
several engines. In the case of a twin-jet air- 
craft the concentration difference ratio distribu- 
tion is given by 



where y and z are respectively the horizontal 
and vertical distances from the aircraft center- 
line, and ECL Is the horizontal separation between 
engines. When the spreading of the exhaust from 
each engine Is small with respect to the distance 
between engines, there Is no significant Interac- 
tion. When the half-value radius approaches the 
engine separation distance, the concentration down- 
stream from each engine Is Increased from what It 
would be alone due to the contribution from the 
other engines. As r^ becomes much larger than 
the engine separation, the concentration downstream 
from each engine Is nearly equal to the number of 
engines times the concentration from each. 


where (c^, - c<»)/(cj - c„) Is the concentration dif- 
ference ratio on the engine centerline for each 
engine. Similar expressions may be formulated for 
other multiple engine aircraft configurations. 

The solution method used In the region where 
atmospheric effects are dominant Is the same as 
the method used In the near wake except that equa- 
tion (10) Is used in place of equation (4) , and the 
atmospheric scale-dependent eddy viscosity Is used 
In place of the turbulent jet eddy viscosity. The 
atmospheric eddy viscosity is formulated using the 
diffusion coefficients given in Reference 3, namely 
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( 12 ) 


where 

2 

K = K = 1 for vertical dispersion 

m z sec 

I 

Kjjj = Ky = 5x10^ for horizontal dispersion 

D = 2x10“^ m^/sec^; atmospheric dispersion rate 

o = for vertical diffusion 

o = y^l’/ZK for horizontal diffusion 

The diffusion coefficients given by Equation (12) 
represent the kinematic eddy viscosity. Thus 
e = pK. 


The dispersion described above Is controlled 
by aircraft and engine effects. The eddy viscosity, 
e. In this region Is decreasing with Increasing 
distance from the aircraft (see Eq. ( 8 )), and at 


The plume width varies continuously from the 
jet dominant region into the atmospheric dominant 
region. In this latter region the plume width ap- 
proaches the solution obtained by assuming scale 
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depeadent dispersion from a point source moving 
with the aircraft. That Is, , 
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for vertical dispersion, and 
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for horizontal dispersion. This solution Is simi- 
lar to that given In Reference 3. In this region 
the horizontal and vertical plume widths are Inde- 
pendent of the aircraft configuration, flight 
speed, the engine power setting (afterburning or 
nonafterburning) , and the number of engines . 


dispersion Is shown on Figure 2. The centerline 
concentration corresponding to the plume size may 
be obtained using Equation (IS) and the Initial 
conditions given In Reference 2. 

In Reference 2, concentrations of CO, NO^, and 
O 3 measured In the wake of an F104 aircraft at dis- 
tances from 2 to 10 kilometers behind the F104 are 
presented and compared with concentrations calcu- 
lated using the near wake dispersion model. Con- 
centrations for this case calculated using the 
present method (with atmospheric dispersion effects 
Included) are Identical to those shown In Refer- 
ence 2 over most of the range of the flight experi- 
ment. Only at the largest distances, -10 km, are 
the calculated concentrations different from those 
given In Reference 2. At that distance, the region 
of atmospheric effects has just begun, and the dis- 
persion Is slightly more rapid than predicted with 
the near -wake model. 


Since In Che far-wake the velocity difference 
Is small compared Co the flight speed, and the 
density of gases In the plume Is nearly equal to 
the ambient density, equations (1) and (10) provide 
the following relation between the centerline con- 
centration difference ratio and the vertical and 
horizontal half -value widths. 


Calculations have been performed for exhaust 
dispersion and dilution for a YF-12 aircraft flying 
at Mach 2.8, 19.8 km altitude, with and without 
afterburning and for a hypothetical four engine 
supersonic transport aircraft of the Boeing SST 
configuration flying at Mach 2.7, 19.8 km altitude 
with maximum afterburning. 


c 


i 



(15) 


For the twin-jet configuration, the concentration 
on the engine centerline Is 
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For an aircraft with (m) engines. In regions where 
the horizontal half-value width Is much larger chan 
Che engine separation distance, the centerline con- 
centration difference ratio and Che horizontal and 
vertical half-value widths are related by 



The horizontal and vertical plume half -value 
widths for these are shown In Figure 3 for the 
YF-12 and Figure 4 for the SST. Also shown is the 
horizontal half -width plus the outboard-engine-to- 
aircraft-centerline distance since this la a meas- 
ure of the overall width of the wake. Atmospheric 
dispersion effects first appear at approximately 
1 minute after aircraft passage. For times on the 
order of an hour or greater, the corresponding 
widths for the YF-12 and the SST are identical. As 
reported by Walton(3) the horizontal diffusivity 
approaches Its asymptotic value In from 5 Co 
10 days; as this occurs the horizontal width growth 
becomes proportional to (t)l /2 indicative of the 
constant diffusivity. Since the asymptotic value 
of the vertical diffusivity is several orders of 
magnitude smaller than Che horizontal value, the 
vertical diffusion Is never really scale-dependent 
and the jet solution, width '(t)l/3, blends 
smoothly into the constant diffusivity solution, 
width -(t)l/ 2 . 


The Independence of the plume size on Che air- 
craft Initial conditions Is an Important result be- 
cause to establish the dilution at any point In the 
far-wake region, only the plume size relations. 
Equations (13) and (14) , and the plume size- 
concentration relation. Equation (17),- are required. 
In the near-wake region, where the plume size Is 
not Independent of the Initial conditions, a dif- 
ferent plume size and/or concentration difference 
curve Is required for every aircraft and flight 
condition of Interest. 

Results and Discussion 


The transition In the plume width from the 
near wake (jet dominant) to the far wake (atmos- 
pheric dominant) Is shown In Figures 1 and 2 for 
Che engines and flight conditions examined In Ref- 
erence 2 . In addition to the curves for the large 
(GE4) and small (J 8 S) engines the solution for the 
plume width for dispersion from a point source 
moving with the aircraft Is given. Vertical dis- 
persion Is shown on Figure 1 whereas horizontal 


The concentration difference ratios as a func- 
tion of time since aircraft passage for the YF-12 
and SST are shown In Figure 5. For times less than 
1 minute the dispersion Is dominated by jet effects 
and the centerline concentration difference ratio 
decrease Is rapid initially, blends Into the asymp- 
totic (t)“2/3 result, and tends toward a t“l/2 
decay as the exhaust jets merge. More rapid dilu- 
tion occurs when the atmospheric dispersion Is dom- 
inant as the eddy diffusivity Is Increasing with 
time. For times from a few minutes to a few days 
after aircraft passage, the concentration differ- 
ence ratio Is proportional to t”^. When the hori- 
zontal diffusivity reaches its asymptotic value, 
the concentration difference ratio Is proportional 
to 1 /t. 


It has been suggested that the turbulent dif- 
fusion In the stratosphere Is not as great as rep- 
resented by the dissipation rate and asymptotic 
diffusivlties used by Walton. ^ ' In Particular Lilly 
(cited In Ref. 4). suggests D = 2x10”^ m^/sec^ and 
Hilst(5) suggests Ky = 10 m^/sec and 
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Kj = 0.1 m^/sec. This value of Ky removes the 
contributions of "meandering," i.e. , turbulence of 
a scale greater than the width of the plume. It is 
not within the scope of the present paper to sug- 
gest values appropriate for stratospheric disper- 
sion, however, it is Interesting to note the sig- 
nificant differences in the predicted dispersion 
obtained by using the above parameters Instead of 
those used in Reference 3. 

In Figure 6 the calculated plume half-value 
widths for the YF-12 at Mach 2.8, 19.8 km altitude, 
with maximum afterburning are shown for both sets 
of stratospheric dispersion parameters. The 
smaller atmospheric dissipation rate and asymptotic 
diffusivitles result in prediction of smaller plume 
dimensions in the region where atmospheric effects 
predominate. Also, the effect of the smaller dis- 
sipation rate is to shift the onset of atmospheric 
effects to larger times after fly-by. 

The centerline concentration difference ratios, 
for the YF-12 flight condition, calculated using 
both sets of stratospheric dispersion parameters 
are shown in Figure 7. Significantly less dilution 
is predicted using the smaller dispersion param- 
eters; for example, at one day after aircraft pas- 
sage the centerline concentration difference ratio 
predicted using D = 2xl0“5 m^/sec^, Ky = 10 m^/sec, 
and Kj = 0.1 m^/sec is 240 times the value pre- 
dicted using D = 2x10"^ m^/sec^, Ky = 5x10^ m^/sec, 
and K^ = 1 m^/sec. 

The concentration of any desired species on 
the plume centerline can be obtained from the 
centerline concentration difference ratio and the 
known (or estimated) ambient concentration and the 
concentration of the species in the engine exhaust. 
Off-centerline concentrations can be obtained also 
by using the half-value width relations in addition 
to the above since the distributions are assumed to 
be Gaussian. 


Summary of Results 

A method is presented for calculating the dis- 
persion of any conserved species in Jet aircraft 
exhaust plumes at cruise conditions. The method is 
appropriate for estimating the dilution of the ex- 
haust from aircraft passage through times on the 
order of weeks thereafter. 

In the near wake of the aircraft, the mixing is 
dependent on the aircraft initial conditions, and 
the flow from each engine is modeled as a round jet 
in a parallel flow. The velocity and concentration 
profiles are assumed to be Gaussian, and the eddy 
viscosity is assumed to be directly proportional to 
the engine airflow plus the absolute value of the 
entrainment mass flow and Inversely proportional to 
the jet width. The concentration decay is rapid 
initially, but approaches the solution for a small 
excess jet as distance from the aircraft Increases. 

Multiple engine aircraft are treated by super- 
imposing the solutions from the several engines. 

At some downstream distance, the eddy viscosity 
calculated for the Jet flow becomes less than the 
eddy viscosity due to scale-dependent diffusion in 
the atmosphere. When this occurs, atmospheric ef- 
fects are assumed to be dominant and the scale- 
dependent eddy viscosity is used. The solution in 


this region approaches that for dispersion from a 
point source moving with the aircraft. Since hori- 
zontal dispersion is known to be more rapid than 
vertical dispersion, a bivariant Gaussian distribu- 
tion for species in the plume is assumed. 


In the far wake region of the flow, the plume 
size is Independent of aircraft configuration and 
engine operating conditions. The dilution at any 
location can be determined from the plume size, 
engine airflow, number of engines, flight speed, 
and altitude. Calculations for exhaust dispersion 
and dilution are presented for a YF-12 and a hypo- 
thetical SST aircraft at high altitude supersonic 
flight conditions. 
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Appendix - Symbols 

centerline concentration - single engine 

concentration on engine centerline for mul- 
tiple engine aircraft 

jet exhaust concentration 
ambient concentration 

2 3 

atmospheric dissipation rate [=] m /sec 

distance between engines for multiple engine 
aircraft 

2 

dispersion coefficient [=] m /sec 

stratospheric dlffusivlty, = K for vertical 
diffusion, = K for horizontal diffusion 

number of engines 

radius 

jet exhaust radius 

half-value radius where (c - c^)/(c - c^) = 

0.5 

plume centerline velocity 
jet exhaust velocity 
streamwise velocity 
flight speed 

time after aircraft passage, t = x/u^ 
streamwise distance 

horizontal distance; (y = 0) = aircraft center- 
line for multiple engines, = engine center- 
line for single engines 

horizontal half-value width 
vertical distance 
vertical half-value width 
eddy viscosity [=] kg/ (msec) 

In 2 

gas density 

plume width standard deviation = r j/->/2A 
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Figure 1. - Exhaust plume vertical half-width for GE4 
and J85 engines at several flight conditions: 

D = ZxlO'^M^/sec^; = 1 Nrlsec; A/B = maximum 
afterburning; MIL= maximum power without after- 
burning. 



Figure 2. - Exhaust plume horizontal half-width for GE4 and 
J85 engines at several flight conditions: D = 2xlO'^M^/sec^; 
Ky= 5x10^ M^/sec; A/B = maximum afterburning; MIL= maxi- 
mum power without afterburning. 




TIME AFTER AIRCRAFT PASSAGE, t, SEC 


Figure 3. - Horizontal and vertical plume widths for YF-12 
aircraft 19.8^ altitude, Mach 2.8: D = 2x10'^ Nrlsec^; 
Ky = 5xl(F M^/sec; = 1 MWsec. 



TIME AFTER AIRCRAFT PASSAGE, t, SEC 

Figure 4. - Horizontal and vertical plume widths for 
SST, 19. 8 kM altitude, Mach 2.7, maximum after- 
burning: D = 2x10'^ M^/sec^; Ky = 5 x 105 M^/sec; 
Kj . 1 M2/sec. ^ 
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Figure 5, - Centerline concentration difference ratio for 
YF-12and SST, 19.8 kM altitude: D = 2x10'^ M^/sec^ 
Ky = 5x10^ M^/sec; = 1 M^/sec. 
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parameters on predicted centerline concentration; flight 
Figure 6. - Effect of assumed stratospheric dispersion condition - YF-12, Mach 2. 8, 19. 8 kM, maximum after- 
parameters on predicted plume widths; flight conditions - burning 

YF-12, Mach 2.8, 19.8 kM. 



